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In the original work by Burnett the pressure tensor and the heat current contain two time 
derivates. Those are commonly replaced by spatial derivatives using the equations to zero order in 
the Knudsen number. The resulting conventional Burnett equations were shown by Bobylev to be 
linearly unstable. In this paper it is shown that the original equations of Burnett have a singularity. 
A hybrid of the original and conventional equations is constructed which is shown to be linearly 
stable. It contains two parameters. For the simplest choice of parameters the hybrid equations have 
no third derivative of the temperature but the inertia term contains second spatial derivatives. For 
stationary flow, when terms Kn 2 Ma 2 can be neglected, the only difference from the conventional 
Burnett equations is the change of coefficients UJ2 —> ^73, ^3 — > 'OJz ■ 
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I. INTRODUCTION 

The practical interest of extending the Navier-Stokes 
equations to smaller length scales comes from the need 
to model re-entrance of space vehicles and in later times 
from nano scale technology. In this parameter region 
there is in general the Boltzmann equation. As the col- 
lision integral is very heavy to calculate numerically, the 
Boltzmann equation is usually applied for fairly simple or 
idealized problems, see Cercignani Q. An other method 
is that of DSMC due to Bird 0. But in the collision 
dominated region also DSMC is very demanding com- 
putationally. Hence special methods are valuable when 
the Knudsen number is too large for the Navier-Stokes 
equations to apply but still small enough to allow for 
an expansion. A fruitful technique in this region is the 
asymptotic method, originating with Hilbert and Grad, 
see Sone 0. 

The Navier-Stokes equations were derived from the 
Boltzmann equation by the Chapman-Enskog method, 
see Chapman & Cowling 4]. They are to first order in 
the mean free path. The corresponding equations to sec- 
ond order in the mean free path were derived by Burnett 
@, see also Q. However, the Burnett equations were 
proven by Bobylev to have a nonphysical instability. 
See also the review by Agarwal et al and the paper by 
Uribe et al. @. 

For the Burnett equations the trivial state of rest is 
thus unstable for perturbations of a wavelength of the 
order of the mean free path and shorter. The Chapman- 
Enskog method is an expansion in the Knudsen number 
Kn = l/L, where I is the mean free path and L a char- 
acteristic length. Wavelengths of the order of the mean 
free path and larger correspond to an effective Knudsen 
number of the order 1 and larger. So there is no contra- 
diction in the fact that the Burnett equations are unsta- 
ble for short wavelengths. The physical content of the 
Burnett equations is for solutions with a characteristic 



length which is large enough. But nevertheless the equa- 
tions should be wellbehaved for all length scales. This is 
necessary mathematically as well as numerically. 

As pointed out by Agarwal et al Burnett's origi- 
nal expression for the viscous pressure tensor contains 
the time derivative of the traceless rate of deformation 
tensor. Chapman & Cowling replace the time derivative 
by spatial derivatives using the equations at the Eulcr 
level. Agarwal et al give the name conventional Burnett 
equations to the resulting equations. There is a similar 
replacement for the time derivative of the temperature 
gradient. This method of replacement using the zero or- 
der equations is part of the procedure in the Chapman- 
Enskog method to obtain the Navier-Stokes equations 
from the Boltzmann equation to first order of the Knud- 
sen number. 

Agarwal et al also raise the question if the instability 
of the conventional Burnett equations is caused by this 
replacement, so that possibly the original Burnett equa- 
tions are linearly stable. In the present paper we show 
that the original Burnett equations have an unphysical 
singularity. In ,9j Jin and Slemrod introduce P, q as new 
independent fields and propose 13 equations, first order 
in time and second order in space. See also the two mas- 
ter theses by Svard flfll and by Stromgren ^lj and the 
paper by Soderholm |l2( based on essentially the same 
idea but with equations which are first order in time as 
well as space. In the paper ^1 these equations are ap- 
plied numerically for nonlinear sound waves. 

In the present paper we introduce a two parameter par- 
tial replacement of the mentioned time derivatives and 
show that the parameters can be chosen to give linear 
stability. We then choose the parameters such that the 
resulting equations are as close to the original and con- 
ventional equations as possible. The resulting equations 
are more similar to the already studied Burnett equations 
than those in and 0] and it seems that it should be 
simpler to reduce a numerical scheme for the equations 
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to one for the Navier-Stokes equations. 



Hence, 



II. THE BURNETT EQUATIONS 

The general equations of balance are 
Dp 



Dt 



Dw 



pV-v = 0, (1) 
-V-P, (2) 
-P : Vv - V • q. (3) 



P- 



Dt 

3fc B DT 
} l>m~Dt 



P is the pressure tensor and q the heat current. We are 
using dyadic notation. 

Let us first write down the original Burnett expression 
for the pressure tensor, see (p = kspT/m) 



P = pl-2^S + n7i — (V-v)S 

P 

+n7 2 ^[^-2(S.(Vv))] 
p Dt 

pT ppT 

+ W ^(VTVT) + w 6 — (S ■ S). 
pi z p 



(4) 



Here. 



S = (Vv) = i[Vv + (Vv) T ]-i(V-v)l. 

1 is the unit tensor, (...) means the symmetric traceless 
part. All other quantities have their usual meaning. The 
original Burnett expression for the heat current is 



q Q = -kV T + 0i-^-(V-v) VT 



(Vv) ■ VT] 



PT 
pT 1 Dt 

9 3 ^!s- Vp + 04 — V • S+6 5 ^—S ■ VT. 
pp ' p pT 



(5) 



In Chapman & Cowling 4] the time derivatives D/Dt 
are replaced by spatial derivatives using the zero or- 
der equations. The resultting expressions are denoted 
Do/Dt. 

To start with we have 

— (VT) = V— - (Vv) • VT. 

r;r ' Dt v ' 

The energy equation to zero order is 
3D T 



2 Dt 



-T(V-v). 



A) 



(VT) 



(G) 



|tV(V-v)-|(V-v)VT - (Vv)-VT. 



Similarly we have 

Dw D r k2 

V l« = T« VV+ ( VV ) < 

< V 7*> = Z* + « W >>- 
The zero order momentum equation is 
D v 1. 



(7) 



TJi 



-Vp. 



Hence, 



^=-<V(-Vp))-(Vv) 2 ). 
Dt p 



(8) 



Explicitly, 



+ ((Vv) 2 ) (9) 

= -^(-VTV/9- -^VpVp + -VVp + VVT). 
pT p p^ p 

Interpreting the time derivatives in |0J and ©) accord- 
ing to © and @ we find the conventional expressions 
for the pressure tensor and heat current. Let us denote 
them P c , q c . The equations of balance CD - © then give 
the conventional Burnett equations. 



III. STABILITY AND REPLACING TIME 
DERIVATIVES BY SPACE DERIVATIVES 



We now make the replacements 



DS DS 

— > (1 — a) Ha- . 

Dt y ' Dt Dt' 



(10) 



D 
Dt 



{VT) ~* (1-/3)^(VT) + /3^(VT), (11) 



where a, [3 are coefficients for which we later shall obtain 
bounds. The choice a — (3 — gives the original Burnett 
equations and a = (3 = 1 the conventional Burnett equa- 
tions. Let us now for simplicity denote nonlinear Burnett 
terms by dots. 

a 2 Dv u 2 
P = pi - 2pS + n7 2 ^(V — ) + tz7 3 -^;(VVT) + ... 
p Dt pi 



This gives 

P = pi - 2pS + m 2 ^(l - a)(V^) 
p Dt 

2 2 

+ -^(w3-atu 2 )(VVT)-an7 2 ^-(VVp) 
pT p z 



(12) 
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In order to study the linear stability of the resulting 
equations we now linearize the equations around a state 
at rest with constant temperature and density. 

p[l + Wa £(l-«)(±A + ±W)]£ 

= -P(l-azu2 — ^A)Vp 

~ f [1 + w (W3 _ atU2 ^ A l VT + ^[Av+iv(V • v)]. 

All the undifferentiated quantities are here taken at the 
background state. For a plane wave with wave number k 
the longitudinal component of the momentum equation 
gives 

dt 



l + n72 g(a-l)§fc 2 mp 

1 + ^( an7 2~ n7 3)|fc 2 fcg 

1 + ^^(^-1)1^ m 



■Viip 



V,,T 



+ 



4/x 



1 + K72 ^(a-I)|fc2 3 P 



Awn 



Putting fc = here we have the linearized Navier- 
Stokes equations. To avoid a singularity in the coeffi- 
cients for a general k we see that it is necessary that 
a > 1. As a = corresponds to the original Burnett 
equations, we see that they have an unphysical singular- 
ity. - For a given k we can interpret the equations as the 
linearized Navier-Stokes equations for an ideal gas but 
with different properties of the gas and the background 
state: 



T 

rh 

1 

TO 



mi + ro2 g(a-l)§fc2 : 

1 l + g(«n7 2 -TI73)ffc 2 

m i + K72 g (a _i)| fc 2 



H = (j, 



This interpretation is possible if the coefficients are pos- 
itive. We see that this requires 

an7 2 — tz73 > 0. 

As < vj2 < ©3 for Maxwell molecules as well as hard 
spheres, a then has to be larger than 1, which excludes 
the original Burnett equations as well as the conventional 
ones. 
Using 

V .S=iAv+iv(V-v) 



we obtain for the heat current 

q=-«VT + , 2 4(l-/^ (VT) 



Dt 



+ 7 [( y )Av + ( y- /3 ? )v(v ' v)] + - 

The linearized energy equation is then 

^ (1+fe ^ (1 _, )A) f 

2m 3 pp dt 

= -p(1 + ^(0 4 -/?02)A)(V-v)+kAT. (13) 

For a plane wave we obtain 



dT 
~dt 

K 

+ - 



T l + ^§W-g 4 )fc 2 
rac « l + 2 I^(/3-l)fc 2 



(V-v) 



A T. 



^l + fcf^/J-l)* 2 
Here we have introduced the specific heat 



Cy 



'ik B 
2m ' 



To avoid singularities it is necessary that (3 > 1. This 
means that the conventional Burnett equations have no 
singularity in the energy equation, but the original Bur- 
nett equations do. We can interpret the energy equation 
as the energy equation for the Navier-Stokes equations of 
an ideal gas with Fourier's expression for the heat cur- 
rent. 



T 

mc v 

k 
c v 



rnc 



• 3 pp 
1 



ci + fcf^Gs-i)^ 



In all we have then given new formal values to the 
five properties of the gas and its background state 
to, c v , (i, k, T. p and v are unchanged. This interpre- 
tation is possible as long as the coefficients are positive. 
As < 64 < 02 for Maxwell molecules as well as hard 
spheres this follows from (i > 1. As the linearized equa- 
tion of continuity is the usual one we conclude that for 
each value k there is an ideal gas such that its linearized 
Navier-Stokes equations coincide with the full set of lin- 
earized hybrid Burnett equations for the same value of 
k. This means that the equations are linearly stable as 
long as 



«>^,/3>l. 

W2 



(14) 



So we have found a two-parameter family of equations 
with linear stability. In this way of thinking, it is easy 
to understand that the conventional Burnett equations 
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are linearly unstable. For large enough k a temperature 
gradient gives rise to a pressure gradient in the opposite 
direction. So if a local region of the order of the mean 
free path or smaller is hotter than its surroundings, the 
pressure in this area will be lower and the gas will flow 
into it, increasing its temperature. But let us stress that 
this is for length scales of the order of the mean free 
path and smaller, so that the equations are not physically 
valid anyhow. As mentioned in the introduction, the 
equations nevertheless need to be well-behaved for such 
perturbations. 

For the transverse components we find 

1 , 21 <9v_l 



P [l + vj2 — {a-l)-k z ] 



PP 



dt 



= /iAv_i_ 



As long as a > 1 this is simply diffusion, with a diffusivity 
depending on k. 



IV. HYBRID BURNETT EQUATIONS 

Let us now make the following choice of the parameters 



a = ,p = 1. 

ZU2 



(15) 



The heat current then is the conventional one, q c . The 
resulting expression for the viscous pressure tensor is now 
denoted P^. 

P h = pi - 2/iS + roi^-fV • v) S 

P 

p 2 DS 
P P 

2 1 rp 

+^3^<— VTVp+-V P Vp-4( Vv ) 2 > 
pT p p* pT 

2 2 

+Wi p^T ( v pV T) + m5 pT^ ( v TV T) 

+TXJQ — (S • S). 

p 

Note that the troublesome third derivative of T is absent. 
Further there is a change of sign in front of DS/Dt as 
compared to the original expression 

Now we have our hybrid Burnett equations. 



Dp 



Dt 



+ pV-v = 0, 

P^r = -V-P fc> 



Dt 
3k b DT 

P 2^Dt = -P^Vv-V-q, 



(16) 
(17) 
(18) 



V. LINEAR STABILITY ANALYSIS 

Let us linearize around a uniform state at rest denoting 
now its temperature Tq and density pq. We write 



T = T (1 + T), p = Po (l + p), d = J— ». 



We introduce the dimensionless variables, where the unit 
of length is of the order of the mean free path 



,p j m _ no m 

X X \ I t m 5 v ' — ^ 



po V k B T 



Po k B T 



In the rest of this section stars and tildes are omitted and 
subscripts denote partial derivatives. We obtain in the 
longitudinal case 

Pt + v x = 
2 

Px ~^-^ZPxxx 

4 2 

+Vt - -^V xx - g( ro 3 - ^2)Vxxt + T x = 

2 4 

% Vx + g(® 4 ~ 02)Vxxx + T t - fT xx = 

/ = (2ro«)/ (3ksp) — b/(3Pr), where Pr is the Prandtl 
number. Looking for solutions exjp(ikx + At) , we find the 
determinant 



A»(l + ^fc 2 ) + A 2 [§(«7 3 - ^ 2 )/fe 4 + (| + f)k 2 } 



+A[(^ - i(0 4 - 2 ))fc 4 + V] + (^ 3 /fc 6 + A 4 )- 

Let us in particular consider the asymptotic case k — » 
oo. We find one mode with 

a _ (ro 3 -w 2 )f 2 
A fc . 

tu 3 

There are also two complex conjugate roots 
A± ps ±i. 



3tx7 3 



2(n7 3 - zu 2 ) 



9zuj + 2(ct 3 - cg 2 )[3ro 3 + 2(6> 2 - 4 )] 

6(n7 3 - TJ7 2 ) 2 / 

For Maxwell molecules, see Q 

2 = 45/8, 4 = 3, n7 2 = 2, w 3 = 3. 

For hard spheres, see 01 

2 = 5.826, 4 = 2.416, tu 2 = 2.029, w 3 = 2.415. 

In both cases zu 3 ,?i7 3 — nj 2 ,02 — 04 are positive. This 
means that the mode corresponding to Aq is damped and 
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FIG. 1: Complex growth factor A for < k < 2. Hard 
spheres. 



nonpropagating. It is the entropy mode. The two modes 
A± are damped. They are propagating sound waves. 

We plot the three roots A in the complex plane in Fig. 
1. Here we use the value the Prandtl number 2/3. This 
is the lowest approximation in terms of Sonine polyno- 
mial expansion for any interatomic potential and is ex- 
perimentally found to be a good approximation, see 0. 
Prandtl numbers in the range between 1/2 to 1 give qual- 
itatively the same plot as Fig. 1. We see in Fig. 1 that 
the real part of A is negativ. 

We have already in the preceding section found that 
the transverse mode is damped and nonpropagating. 



VI. DISCUSSION 

Let us pull out the acceleration from the viscous pres- 
sure tensor P/j using ) 



P h = {w 2 - W3 )^((V— '-)) + P h 
p Dt 



Here 



P h =pl-2/xS + TOi — (V-v)S 

P 

_( Wa - roa ) — <(Vv) 2 ) 
P 

-2 ro2 ^(S-(Vv))-tn 3 4(VVp) 
P P 

2 1 rp 

pT p p* pT 

2 2 

+^4 (V PV T)+m 5 ^ (V TV T) 
ppl pi z 

+ru e — (S ■ S). 
P 



The momentum equation can now be written 



^-V-[(-3-^((V^)>] 



VP 



h- 



Let us write a for the acceleration and study the terms 
containing the acceleration 

u 2 

pa- V- [(tj7 S - zu 2 ) — (Va)] 
P 

We multiply by a and integrate over the region of flow 

/{a • [pa - V • ((tu 3 - w 2 ) — <Va))]}dF 
J P 

= J[pa 2 + (vj 3 -vj 2 )^-(\/a):{Va)}dV 

- f(w 3 - w 2 )— a ■ (Va) • dS 
J P 

If there are no boundaries and the flow vanishes at in- 
finity, the surface integral vanishes. We conclude that 
the operator acting on Dv/Dt is then positive and the 
momentum equation can be solved for the acceleration. 
The situation when there are boundaries requires a closer 
examination. 

In order to see more clearly the structure of the hy- 
brid Burnett equations, we replace the nonlinear terms 
in Kn 2 by dots. But first we split up in longitudinal and 
transverse parts according to 

iA + iw = ^VV + i(A-VV). 

The hybrid Burnett equations can then be written 

Dp 



Dt 



+ p(V-v) = 0, (19) 



p{l - (H7 3 - ^)i-[\{Al - VV) + |VV} • ^ (20) 
pp 2 3 Dt 



u 2 2 

-Vp + 2V • (p,S) + w^-AVp - 
p z 3 



-^p— = -p(V ■ v) + 2pS ■ S + V • (k V T) (21) 
2m Dt 

-(04-0 2 )-~A(V-v) + ... 
p 3 

Let us also write down the momentum equation of the 
conventional Burnett equations There is no need to write 
down the equations of continuity and energy as they are 
exactly the same as for the hybrid equations. 



-Vp + 2V-(pS)-(^ 3 -^ 2 )~AVT (22) 

pi 3 



p 2 2 

+W2^-AV/5- 

p z 3 



The difference is that the hybrid Burnett equations have 
a more complicated inertia term but no term A VT in the 
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momentum equation. The coefficient in front of the term 
AVp is w 3 — ZU2 in the conventional Burnett equations 
but 073 in the hybrid Burnett equations. 

Let us also note that the equations without the Kn 2 
nonlinear terms are valid when Kn 2 Ma 2 can be neglected 
and the relative temperature and density variations are 
assumed to be of the order of Ma. 

Let us take a closer look at the inertia term in the 
hybrid Burnett equations, see (|T7jl . We consider a Fourier 
component of the acceleration Dv/Dt proportional to 
exp(ik • r). Then 

p{[l - (w 3 - ro2 )£-[~(Al - VV) + |VV]}^ 

pp 2 6 JJt 

2, N u 2 fc 2 1 kk lJ Dv 

We have effectively a transverse inertia in the first term 
and a longitudinal inertia in the second term. As 
ru 3 — m-i > 0, both of them are positive. We noted earlier 
that the sign of the Dv/Dt term of the hybrid Burnett 
expression for P is the opposite of that of the original 
Burnett equations. As a result the original Burnett equa- 
tions have a singularity for certain wave numbers where 
the inertia vanishes. 

Let us also consider the low Ma stationary case. The 
hybrid Burnett equations are then 

V • Gov) = 0, 



p(v • V)« = -Vp + 2V • (pS) + m 3 t--AVp, 

p z 3 

— — p(v-VT) = -p(V • v)+2^S : S 
2m 

+V-( K VT) - (0 4 -0 2 )£LHa(V-v). 

9 3 

The conventional Burnett momentum equation is 

j o(v-V)v = -Vp + 2V-(pS) (23) 

u 2 2 u 2 2 

-(w 3 -tn 2 )4--AVT + ti7 2 ^-AVp, 
pi 3 3 

The only difference between the hybrid Burnett equa- 
tions and the conventional Burnett equations is the 
change of coefficients 

W 2 — > 073, TU 3 — > V3 3 . (24) 

This also means that the third derivative of the temper- 
ature is absent in the hybrid Burnett equations. 
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